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We propose a general construction principle which allows 
to include an infinite number of resonance states into a scat- 
tering matrix of hyperbolic type. As a concrete realization 
of this mechanism we provide new S-matrices generalizing a 
class of hyperbolic ones, which are related to a pair of simple 
Lie algebras, to the elliptic case. For specific choices of the al- 
gebras we propose elliptic generalizations of affine Toda field 
theories and the homogeneous sine-Gordon models. For the 
generalization of the sinh-Gordon model we compute explic- 
itly renormalization group scaling functions by means of the c- 
theorem and the thermodynamic Bethe ansatz. In particular 
we identify the Virasoro central charges of the corresponding 
ultraviolet conformal field theories. 

FAGS numbers: 11.55.Ds, ll.lO.Hi, ll.lO.Kk, 05.70.Jk 



I. INTRODUCTION 

Treating quantum field theories in 1-1-1 dimensions as a 
test laboratory for realistic theories in higher dimensions, 
this paper is concerned with the general question of how 
to enlarge a given finite particle spectrum of a theory to 
an infinite one. 

In general the bootstrap ||l|, which is the construction 
principle for the scattering matrix, is assumed to close 
after a finite number of steps, which means it involves 
a finite number of particles. However, from a physical 
as well as from a mathematical point of view, it ap- 
pears to be natural to extend the construction in such 
a way that it would involve an infinite number of parti- 
cles. The physical motivation for this are string theories, 
which admit an infinite particle spectrum. Mathemati- 
cally the infinite bootstrap would be an analogy to infi- 
nite dimensional groups, in the sense that two entries of 
the S-matrix are combined into a third, which is again a 
member of the same infinite set. It appears to us that 
it is impossible to construct an infinite bootstrap system 
involving asymptotic states and find the mathematical 
analogue to infinite groups in this sense (see also foot- 
note 2 and the paragraph after figure 1). However, it 
is possible to introduce an infinite number of unstable 
particles into the spectrum. Scattering matrices which 
would allow such type of interpretation have occurred in 
the literature although only in the latter paper a 

reference to unstable particles has been made. In 
these matrices were found to be expressible in terms of 
elliptic functions, a feature very common in the context of 
lattice models, e.g. Q. The main purpose of this paper is 
to suggest a general construction principle for such type 



of S-matrices starting from some known theory with a 
finite particle spectrum of a special, albeit quite generic, 
form. As particular examples we provide elliptic gener- 
alizations of scattering matrices related to a pair of Lie 
algebras |^ , which contain the affine Toda S-matrices |^ 
and homogeneous sine-Gordon S-matrices ^ for partic- 
ular choices of the algebras. 

Our paper is organized as follows: In section II we 
provide a general principle for the construction of scat- 
tering matrices which involve an infinite number of un- 
stable particles and present some explicit examples. In 
section HI we construct renormalization group (RG) scal- 
ing functions by means of the c-theorem and the thermo- 
dynamic Bethe ansatz (TEA) for the generalization of 
the sinh-Gordon model. Our conclusions and an outlook 
towards open problems are stated in section IV. 



II. CONSTRUCTION PRINCIPLE 

Let us consider the huge class of two-particle S- 
matrices, which describe the scattering between particles 
of type a and & as a function of the rapidity difference 9, 
of the general forin|^ 



(1) 



Here S^^{9) denotes the so-called minimal S-matrix 
which satisfies the consistency relations Q, namely uni- 
tarity, crossing and the fusing bootstrap equations and 
possibly possess poles on the imaginary axis in the sheet 
< lm6 < TT, which is physical for asymptotic states. 
The CDD-factor §, referred to as S^^^{e), also sat- 
isfies these equations, but has its poles in the sheet 
— TT < Im6' < 0, which is the "physical one" for resonance 
states. S™^{9) might depend on additional constants 
like the effective coupling constant or a resonance param- 
eter. A simple prescription to introduce now an infinite 
number of resonance poles is to replace the CDD-factor 
in (0) by 



N 



(2) 



^Exceptions to this factorization are for instance the scatter- 
ing matrices of affine Toda field theories related to non-simply 
laced Lie algebras, which was first noted in llOl. 
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where lo is taken to be real. By construction the new 
S-matrix, Sabi9,N) = S^^"" (6) S^^^ {9 , N) satisfies the 
bootstrap consistency equations and possible poles in the 
sheet — TT < Im 6* < have now been duplicated 2N times 
within this sheet, such that they admit an interpretation 
as unstable particles. Therefore, when N ^ oo we have 
an infinite number of resonance poles. Since, as a conse- 
quence of crossing and unitarity, the S-matrix is known 
to be a 27ri-periodic function, a property shared indi- 
vidually by S^^^{9,N), we expect to recover a double 
periodic function in the limit N oo 



N 



hm S^ri9,N)=\im S^^i 



N- 



+ fi2Tri + iyuj,N) (3) 



for /X, ^ g Z. At this stage it is not clear whether the 
prescription (||) is meaningful at all, in the sense that it 
leads to meaningful quantum field theories, and in par- 
ticular one has to be concerned about the convergence 
of the infinite product in (||). Since limjv^oo S^b'^{9,N) 
is a double periodic function we expect that it is some- 
how related to elliptic functions (see e.g. for their 
properties). Let us therefore now look concretely at the 
building blocks which can be used to make up the entire 
scattering matrix in the non-elliptic case, when backscat- 
tering is absent. In that case the S-matrices are diagonal 
and known to be of the general form 



Sabi9) ^ Y[{X}\ 



ntaiih{0 ~ inx + a)/2 



xeA 



with X G Q and cr G R. A specific theory is then char- 
acterized by the finite set j^. This means, if we demon- 
strate that the prescription is meaningful for each 
individual building block {xjg as defined in (Q), in par- 
ticular we need to demonstrate the convergence of the 
infinite product, we have established that it is sensible 
for the entire scattering matrix. For this purpose we 
note the identity 



{^6 



n 



sc 0_ dn ( 
"dn^ 



sc ( 



(5) 



Here we abbreviated 9± = {9 ± inx + a)iKi/TT and used 
the Jacobian elliptic functions in the standard notation 
pq(z) with p,q G {s,c,d,n} (see e.g. 0). The quarter 
periods depending on the parameter £ e [0, 1] are 
defined in the usual way through the complete elliptic 
integral 



7r/2 



{I -tsin^ 9)-^/'^ d9. 



(6) 



The period of {x}gf, is chosen to be w = nK^^i^i)/ Ki. 
The last identity in (||) is easily derived from the infinite 
product representations of the elliptic functions which 
can be found in various places as for instance in [ pl| 

TTo: ^ l-2g^"cos(7ra;/j^,) + g4" 
sc a; - /c tan ^ ^ ^ ^^^^ cos{nx/Ki) + g^n ' 



dnx — 



1 + 2g2"-i cosiirx/ Ki) + g4»-2 



n=l 



- 292n-i cos{7rx/Ki) + q 



4n-2 ' 



(8) 



with k — and q — exp(— w). Recalling the well 

known limits lim^^o ^£ — '^/'^ and lim^^o = oo 
we obtain 



lhn{a;}^^, - {x}? 



(9) 



This means in the limit i ~* the elliptic S-matrix Sab{9) 
collapses to the hyperbolic one, that is Sab{9)- Notice 
that due to the general identity pr(a;)/ qr(a;) — pq(x), we 
could also write (^) in terms of various other combina- 
tions of elliptic functions. For instance replacing sc by 
sn / cn is probably most intuitive, since it allows an alter- 
native prescription to (^ for the construction of elliptic 
scattering matrices: Replace sinh sn, cosh — > cn and 
correct the consistency equations by a factor dn, which 
reduces always to 1 in the hyperbolic limit, in such a way 
that no resonance poles are left inside the physical sheet. 
Defining now the function 

0f,.y{x, cr, uj) := 27ri(t/ + x/2) + 2ujfi - a (10) 

the singularities of {x}gg are easily identified as 



- x,a,uj), (11) 



poles: 9^,y{-x,a,u}), 6^y{x - 1, a,u;) . (12) 

Note that when taking < a; < 1 the poles are situated 
in the non-physical sheet. 

This brings us to the question of how to interpret these 
poles and how can we characterize the physical properties 
of the related particles? Considering the S-matrix Sab{9), 
which describes the scattering of two particles of type 
a and b with masses uia and mb, we assume that there 
is a resonance pole situated a.t 6^ — a ~ ia. According 
to the Breit-Wigner formula jl7| (see also e.g. [|l8|) the 
mass Mc and the decay width of an unstable particle 
of type c can be conveniently expressed as 



2Mi 



+ Fi/2 = Vt^+I^-^, (13) 



where 



The fact that a; £ Q together with the bootstrap leads 
to a finite set A and therefore a finite number of asymptotic 
states. Taking instead a; G R could possibly lead to an infinite 
number, but a consistent closure of the bootstrap is not known 
up to now. 



7 = -h TO^ -|- 2maTOh cosh a cos cr, 
7 = 2mamb sinh | cr | sin ct . 



(14) 
(15) 



We keep here in mind that this description, although 
frequently used e.g. l9p,0], is not entirely rigorous and 
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requires additional investigation. This caution is based 



on various facts. First, the relations (13) are simply de- 
rived by carrying over a prescription from usual quantum 
mechanics to quantum field theory, i.e. complexifying 
the mass of a stable particle. Second, solving the Breit- 
Wigner formula for the quantities A/g, Fg and treating 
them literally as mass and decay width is somewhat prob- 
lematic since this is in conflict with Heisenberg's uncer- 
tainty principle, because apparently we know simultane- 
ously the energy and the time. Third, the Breit-Wigner 
relations presume an exponential decay in momentum 
space, which is in fact incompatible with the general 
principles of quantum field theory and therefore might 
possibly be a problem in this context Nonetheless, 
we employ these quantities and try to find evidence to 
support that they are indeed meaningful. When taking 
Mc to be the mass of the unstable particle there should 
be a threshold for energetic reasons of the type 



(16) 



with the consequence that the decay width is bounded 
by 



F| > 8mami,{l — cosher cos cr) 



(17) 



So far evidence for these thresholds has not been found 
in the literature. One reason for this is that the unstable 
particles enter the bootstrap principle in a more passive 
way than the stable particles, whose properties are di- 
rectly used in the construction procedure. Hence one 
expects that signs for these thresholds will emerge in a 
more indirect way. 

A summary of our statements about the pole structure 
of S{0) is depicted in figure 1. 




Figure 1: The poles of the blocks {x}g ^ are the crosses in the 
sheet — TT < Im < 0. The crosses on the positive part of the 
imaginary axis are associated, as usual, with stable particles. 
For equal masses of the stable particles the threshold ( |l6| ) is 
at — arccosh[(3 — cos7ra;)/(l -I- costtx)]. 



Since the poles inside the sheet < Im0 < tt are as- 
sociated to S™^'^{9), it is also obvious from figure 1 why 
the prescription (||) may not be employed for this part 
of the S-matrix, since it would lead to a pole structure 
which is, according to (O), non physical for and Fg. 



A. Examples 

It is clear by construction that our prescription in- 
cludes all affine Toda field theories related to simply laced 
Lie algebras, since they all factorize as (0) and may be 
represented in the form (^), see e.g. Taking the res- 
onance parameter a to be zero and the set A = {t\ for 
< i < 1, we recover as a special case the elliptic version 
of the sinh-Gordon model proposed in . Reintroducing 
CT, its scattering matrix reads 



s{e)= n 



n— — oo 



t(m\\{9 — inx + nuj + cr) /2 
tanh(0 -|- inx + noj + a) /2 



(18) 



According to (g^ the masses and decay width of the 
unstable particles are 



iV2 cosh 



F;^;'^" =m2\/2sinh 



y 



y 



-x,x-l (19) 
-x,x-l (20) 



where m denotes the mass of the stable particle. The 
thresholds (|6h and (O) translate in this case into 



„ ■ 2 tt(2x + 1±1) 

, 3 =F COS TTCc _^ , sm — j - 

cosh(nw + a) > — , 1 > 4m- 



± cos TTX 



cos 



7r(2a: + l±l) 



(21) 



As a further example we consider the elliptic general- 
ization of the Ai|v4jv_i-theory (= 5't/(A^)2-homogeneous 
sine-Gordon model). The two-particle S-matrix describ- 
ing the scattering of two stable particles of type a and 6, 
with l<a, 6<A^— 1, related to the non-elliptic version 
of this model was proposed in . In our notation it may 
be written as 



SabiO,aab) = (-1)'"" 



{1/2}- 



(22) 



Here / denotes the incidence matrix of the SU{N)- 
Dynkin diagram, the resonance parameters have the 
property aab = —crba and Ca = ±1 depending on whether 
a is even or odd. According to our prescription outlined 
in the previous paragraph, the elliptic generalization of 
©is 



Sabi9,aab,i) = (-1)'"" Ca./{l/2}^; 



(23) 



Note that despite the appearance of the square root, S 
as well as S are still meromorphic functions in 9. 
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III. RG-SCALING FUNCTIONS 

Having established that our prescription leads to sen- 
sible solutions of the bootstrap consistency equations, we 
would also like to know what kind of quantum field the- 
ories these scattering matrices correspond to. Up to now 
all known solutions to the on-shell consistency equations 
have led to sensible QFT's, albeit a rigorous proof which 
would establish that indeed all solutions are well-defined 
local QFT's is still an outstanding issue. Some crucial 
characteristics of the theory are contained in the renor- 
malization group scaling functions, which we now want 
to determine. In particular, we want to identify in the 
extreme ultraviolet limit the Virasoro central charges of 
the corresponding conformal field theories. 



lim c'^-* 



(r) 



de 



cosh" e 



(25) 



It is here crucial to note that besides the formulation 
of S in terms of elliptic functions for iV — > ck), it can 
also be expressed equivalently in terms of the usual sinh- 
Gordon S-matrix (^. When trying to solve now the form 
factor consistency equations ||20| , we can exploit this ob- 
servation. Since for the model at hand there is neither 
a kinematic nor a bound state pole in FfiO), the only 
equations to be solved are Watson's equations. The two 
particle form factor is then easily obtained to be 



N 



(26) 



A. The c-theorem 

We carry out this task by evaluating the c-theorem 
in the version presented in [|l3[ 



oo 



ddi . . . dOn ^rE 

e 



n—1 fj,i ...fj,n_^ 

F«l^-^"(0i,...,0„) 



n!(27r) 

2 (6 + 6rE + 3r^E^ 



(24) 



2E^ 



The sum of the on-shell energies is here denoted by 
E = X]r=i ^fj-i cosh^i, with TO^. being the masses of the 
theory and the correlation function for the trace of the 
energy momentum tensor Q has been expanded in terms 
of n-particle form factors (6*1, . . . , 6'„) (see ||2^] 

for general properties and jl^ for explicit sinh-Gordon 
formulae). We normalized O and m^. by an overall 
mass scale, such that E as well as the renormalization 
group parameter r become dimensionless. In particular 
limj._,o c(r) is the ultraviolet Virasoro central charge. 

Let us now start with the evaluation of c(r) as defined 
in (|2^) for the elliptic version of the sinh-Gordon model. 
As the input for this we need to know the n-particle form 
factors. Since so far it is not known how to compute the 
sum in n analytically, we have to resort to a numerical 
treatment and it is clear that we have to terminate the 
series at a certain value of n. Fortunately, it was ob- 
served explicitly in [|l6|, that in fact the expression for 
n = 2 is already very close to the exact answer for the 
sinh-Gordon model. We assume here that the conver- 
gence behaviour is still true when we generalize the scat- 
tering matrix to (|lj) . Note that in general one has to be 
careful with this approximation, since the higher parti- 
cle contributions are crucial in some models in order to 



obtain a good approximation to c(r) ]21|, 13| , p2| . In the 
two-particle approximation, indicated by the superscript, 
one can perform one of the integrations analytically and 
(p4) acquires the simple form 



where -F'inin(^) is the minimal form factor of the sinh- 
Gordon model obtained in IM 



^'min(6') = exp 



L 



sm 



h(i^)sinh(f )sinh(|) 



sinh(t) 



(27) 



Using the infinite product representation for i^min(^) 0] 
the solution ( p6| ) for N ^ oo coincides with the equation 
(5.3) in Proceeding now to the evaluation of (P^), we 
require {^2* {0 , N)\'^ , whose characteristics are captured 
in figure 2. 
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Figure 2: Absolute value squared of the two particle form fac- 
tors f{8,N) = \Ff{e,N)/27r\'^ as functions of the rapidity 
for different values of N for u) = 1.3 and x — 0.1. 

We observe that for a certain value of TV the function 
starts to converge, which is of course important from a 
technical point of view when we want to compute the 
limiting case N oo. The other observation we make 
in figure 2, see also figure 5 in [H, is that f{0,N = 0) 
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always has a distinct maximum, which we refer to as 9m- 
From ( ^tI) follows that it is determined by the solution 
of 



4a 



cosh 9m sin ttx + cosh 2xt: coth — 

TT 2 

Mm. 

2 



cosh 



sinh 



0„ 



2{2x — 1) cos nx sinh 9„ 



(28) 



Solving this equation for various values of x, we find that 
6m is always slightly greater than the smallest threshold 
bound obtained from For instance for a; = 0.1 we 

obtain 9m — 1-439 and nu; + a > 0.315 and for x — 0.5 
we have 9m — 2.040 and nu + a > 1.763. We interpret 
this as an indication that the form factors "know" about 
the thresholds (|2l|). We support this now by considering 
limAT^oo f(^, N) for various values of uj. 




Figure 3: Absolute value squared of minimal form factors 
g{d,Lu) — limAr_>oo \Fmini9, N)\'^ as functions of the rapidity 
for different values of lo and x — 0.1. 



h(e,(o) 















co=0.2 








O)=0.25 

£0=0.5 

co= 1.0 






m 














Figure 4: Integrand of equation (|25[), that is h[9,uj) = 
limjv^oo |^inin(^, cosh^ as a function of the rapidity 

for different values of lu and x — 0.1. 



We observe that in the region in which the factor 
l/cosh^(6»), emerging in (|2^), is still non vanishing the 
integrals limjv^oo / d9\Fjnin{9, N)]"^ are decreasing func- 
tions of LU. This behaviour is changed once we take 
UJ < 9m &s we can explicitly extract from figure 4. 

Naturally this features are also reflected in the scaling 
functions. Presuming that for each value of N we have a 
consistent theory, we would like to know which ultravio- 
let central charges these models possess and in addition 
we want to identify a value of N for which the related 
model constitutes reasonably good approximation for the 
elliptic models. That such an identification is possible is 
exhibited in figure 5. In addition we observe, that for 
fixed LO and x the scaling function is a monotonically in- 
creasing when N is varied. 
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Figure 5: Ultraviolet Virasoro central charge c^^-* as a function 
of iV. 

Focussing now on the elliptic case, that is we select a 
large enough N such that this case is well approximated, 
we compute the scaling function in dependence of lu for 
various values of r. Our results are depicted in figure 6. 




Figure 6: Ultraviolet Virasoro central charge c^^' as a function 

of UJ. 
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In the extreme limits we obtain limt^^o c(r, w) = 
and for large uj we recover the values of the sinh- 
Gordon model, lirn^^oo c(r, w) = csQ{r). The latter 
limit follows from (|9|) and is in addition compatible with 
limg-^oo -Fniin(^) — 1- Thc values for the extremal points 
were already quoted in however, we also observe 
that the function is not monotonically increasing between 
these points as claimed in there. In fact, in the physi- 
cal region, that is for values of w > 0.315, the function is 
monotonically decreasing and does not take on values be- 
tween and 1. Remarkably, the threshold is quite clearly 
exhibited by a drastic change in the behaviour of c, that 
is the onset of a small plateau as is visible in figure 6. 
We performed the same computation for different values 
of X and observed that this onset moves in the direction 
predicted by equation (|2l|). 



B. The thermodynamic Bethe ansatz 



Using alternatively the representation of the S-matrix 
( |l8| ) in terms of elliptic functions, we compute the kernel 
directly to 



lim ipN{9) = — V 



k=- 



dc Ok 
sn 9k 



snt/fc 



(33) 



With these expression we carry out our numerical anal- 
ysis, that is we solve iteratively the equation ( p9[ ) and 
evaluate (|30| ) thereafter. The results of this investiga- 
tions are presented in figure 7. 

Unfortunately for very small values of uj and r our nu- 
merical iteration procedure does not converge reliably. 
However, we will be content at this stage with the data 
obtained so far, since they already support qualitative 
our c-theorem analysis. They confirm that above thresh- 
old the scaling function is monotonically decreasing as a 
function of uj and also that values greater than 1 may be 
reached, even for finite values of r. 



Let us now compare the results of the previous section 
with an alternative method, namely the thermodynamic 
Bcthe ansatz For this we first have to solve the 

TBA-equations 



rrhi cosh0 + ln(l 



E 

3 



ip,j*L,{0) (29) 



for the function Li{9). The information of the scat- 
tering matrix is captured in the kernel ipij (0) — 
—id\aSij{9)/d9 of thc rapidity convolution, which is de- 
noted as usual by / * g\9) := J d9' /2t: f [9 - 9')g{9'). 
The dimensionless parameter r — miT~^ is the inverse 
temperature T times the overall mass scale of the light- 
est particle mi. Also all masses have been normalized 
in this way, i.e. rhi — mi/mi. Having determined the 
Li(0)-functions, we may compute the scaling function by 
means of 



c'(r) 



3r 



,2 ^ 



d9 cosh9L,i9) 



(30) 
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Figure 7: TBA scaling function. 



IV. CONCLUSIONS 



Once again linv^o c'{r) is the ultraviolet Virasoro central 
charge. We would like to recall here that the scaling 
functions c(r) and c'(r) are not identical, but contain 
qualitatively the same information in thc RG sense. 

In order to carry out this analysis we need to know 
in (29) the kernel (p{9) as input. For the model under 
consideration we can exploit the factorization property 
(p^) for a finite product and trivially obtain 



N 



(9) = ^80(9 + nuj + a) , 



n=-N 

where fsG{9) is the sinh-Gordon kernel, e.g. 

4sin(7ra;) cosh6' 



Vsg{9) = 



cosh(20) — cos(27ra;) 



(31) 



(32) 



Starting from a given scattering matrix of hyperbolic 
type, we have demonstrated that it is possible to include 
consistently an arbitrary number of unstable particles 
into the spectrum of the theory. In particular when this 
number becomes infinite the S-matrix may be expressed 
in terms of elliptic functions. 

For the generalization of the sinh-Gordon model we 
computed RG scaling functions. Within these analysis we 
found clear evidence for the thresholds which constrain 
the masses of the unstable particles. Above threshold, 
the values the ultraviolet Virasoro central charges may 
take are between 1 and 2 (possibly slightly greater than 
2) and not between and 1 as suggested in The 
theories are consistent for each finite value of N. For 
fixed resonance parameters uj and a the scaling functions 
are non-decreasing for increasing N. 
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Concerning the investigation of the c-theorem, it would 
be desirable to refine the analysis. In particular one 
should include higher n-particle form factors into the 
expansion. For the elliptic version some of them were 
already presented in but in general it remains a 
challenge to find closed expressions for arbitrary particle 
numbers. At present the TBA analysis is the least con- 
clusive exploration and deserves further consideration in 
future. In particular the regions of uj and r, which were 
not accessible to us, should be explored and might pos- 
sibly lead to a further more concrete indication of the 
thresholds also in this context. In regard to this, it will 
be useful develop existence criteria for the solution of the 
TBA equations analogue to the one derived in Q . The 
one presented in there can not be taken over directly, 
since it makes use of the fact that J d9\ip{9)\ equals 2tt, 
whereas for the model investigated here this is 2TrN. It 
would be desirable to develop analytic approximations 
for the TBA solutions in the extreme ultraviolet limit, 
i.e. r = 0, similar to the ones already existing for theo- 
ries with different characteristic features. 
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